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ABSTRACT 

We argue that eclipses of radio emission from the miUisecond pulsar A in the double pulsar system 
J0737-3039 are due to synchrotron absorption by plasma in the closed field line region of the magneto- 
sphere of its normal pulsar companion B. Based on a plausible geometric model, A's radio beam only 
illuminates B's magnetosphere for about 10 minutes surrounding the time of eclipse. During this time 
it heats particles at r > 10^ cm to relativistic energies and enables extra plasma, beyond that needed 
to maintain the corotation electric field, to be trapped by magnetic mirroring. An enhancement of 
the plasma density by a factor ~ 10^ is required to match the duration and optical depth of the 
observed eclipses. The extra plasma might be supplied by a source near B through pair creation 
by energetic photons produced in B's outer gap. Relativistic pairs cool by synchrotron radiation close 
to where they are born. Re-excitation of their gyrational motions by cyclotron absorption of A's radio 
beam can result in their becoming trapped between conjugate mirror points in B's magnetosphere. 
Because the trapping efficiency decreases with increasing optical depth, the plasma density enhance- 
ment saturates even under steady state illumination. The result is an eclipse with finite, frequency 
dependent, optical depth. After illumination by A's radio beam ceases, the trapped particles cool and 
are lost. The entire cycle repeats every orbital period. We speculate that the asymmetries between 
eclipse ingress and egress result in part from the magnetosphere's evolution toward a steady state 
when illuminated by A's radio beam. We predict that A's linear polarization will vary with both 
eclipse phase and B's rotational phase. 

Subject headings: pulsars: general — pulsars: individual (J0737-3039A, J0737-3039B) — stars: neu- 
tron — radiation mechanisms: non-thermal — plasmas 



1. INTRODUCTION. 

The recent discovery of the binary pulsar J0737-3039 
— a millisecond pulsar (pulsar A with a period Pa — 23 
ms) and a normal pulsar (pulsar B with a period Pb = 
2.8 s) in a tight 2.4 hrs orbit (Burgay et al. 2003) — 
not only provides us with unprecedented tests of general 
relativity (Lyne et al. 2004) but also reveals a variety of 
magnetospheric phenomena. Among the latter are varia- 
tions of pulsar B's radio emission correlated with binary 
orbital phase (Lyne et al. 2004; Ransom et al. 2004) and 
modulated at the spin frequency of pulsar A (McLaugh- 
lin et al. 2004a), and periodic eclipses of pulsar A when 
it passes behind pulsar B (Lyne et al. 2004; Kaspi et 
al. 2004; McLaughlin et al. 2004b). It is the latter phe- 
nomenon that concerns us in this paper. 

Detailed observations by the Green Bank Telescope 
(Kaspi et al. 2004) established the following frequency- 
averaged properties of pulsar A's eclipses: the eclipse 
duration is about 27 s, which for a relative transverse ve- 
locity of 680 km s~^ translates into a size of 18, 000 km. 
Eclipses are significantly asymmetric with ingress taking 
3 — 4 times longer than egress; pulsar A's radio beam 
is extinguished more strongly post conjunction, consis- 
tent with zero fiux, than before conjunction when some 
flux leaks through. Analysis of the same data at higher 
time resolution by McLaughlin ct al. (2004b) uncovered 
effects of B's rotational phase on the frequency depen- 
dence and shape of the eclipse during ingress; its shape 
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during egress is remarkably independent of either B's ro- 
tational phase or radio frequency. 

The spin-down luminosity of pulsar A exceeds that of 
pulsar B by about 3600, so it is plausible that B's magne- 
tosphere is compressed by a relativistic wind from A. Cal- 
culations by Lyutikov (2004) demonstrate that the ram 
pressure of A's wind can be balanced by B's magnetic 
field pressure at a stand-off distance rso ~ (3.5 — 6) x 10^ 
cm which is within B's light cylinder radius of t^ b = 
c/VIb ~ 1-3 X 10^" cm. A crucial point is that the size 
of the eclipsing region is considerably smaller than even 
the compressed size of the B's magnetosphere. To fully 
quantify the geometry of the eclipse, the inclination of 
the system must be accurately known. Measurements of 
the Shapiro delay established that the inclination is very 
high, i = 87° ± 3° (Lyne et al. 2004). A refined estimate 
by Coles et al. (2004) based on correlation of interstellar 
scintillations of both pulsars yields i — 90.26° ± 0.13°, 
which corresponds to a minimum distance of the A's ra- 
dio beam with respect to B's position projected on the 
plane of the sky of only 4000 ± 2000 km. This refinement 
implies that the extinction during eclipse arises inside 
pulsar B's magnetosphere and that the radial extent of 
the eclipsing region is about 10, 000 km. 

Our goal is to evaluate the absorption of the radio 
beam of pulsar A as it passes through the magnetosphere 
of pulsar B. We show in §21 that resonant cyclotron ab- 
sorption in the charged-separated, closed field line region 
would provide only a small optical depth. However, we 
demonstrate in fj^that absorption of radiation from pul- 
sar A heats charged particles in B's magnetosphere to 
relativistic energies. In ^we describe how this enables 
the accumulation of additional, charge neutral plasma 
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in B's magnetosphere. As a result, the radio emission of 
pulsar A can be significantly extinguished by synchrotron 
absorption. We devote J7| to a discussion of the ramifi- 
cations of our model. 

2. RESONANT CYCLOTRON ABSORPTION. 

The region of closed field lines in the conventional pul- 
sar magnetosphere model contains a corotating, charge- 
separated plasma with number density 



I^B •B(r)| 
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(1) 



where v is the corotation velocity (Goldreich & Julian 
1969). We are primarily interested in the region of B's 
magnetosphere where v c and the magnetic field is ap- 
proximately dipolar. For our simplified model it suffices 
to ignore the angular dependence of the field and set 
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where i?^ is the neutron star radius and Bq is its surface 
magnetic field. 

The particle number density n can be higher than noj 
because the addition of neutral plasma does not affect 
the net charge density. We characterize this increase by 
the parameter A > 1 defined such that 



n{r) « A(r) 
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In what follows, we assume that the magnetospheric par- 
ticles are electrons and positrons. 

Resonant cyclotron absorption is the dominant source 
of extinction for radio waves passing through a conven- 
tional pulsar magnetosphere. Cyclotron resonance oc- 
curs where 

u;KiujB = , (4) 

TTleC 

with cross-section (Canuto, Lodenquai & Ruderman 
1971; Daugherty & Ventura 1978) 



Here (tt is the Thompson cross-section and 
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the natural line width. Strictly speaking, this cross sec- 
tion applies to photons in particular modes which differ 
for electrons and positrons, but we ignore this detail here. 

The peak cross section at resonance reaches amax ~ 
67r(c/ajB)^, or roughly the square of the photon wave- 
length. Nevertheless, the optical depth is modest because 
an incident photon only resonates with cold electrons in 
a narrow radius range. From equations we find 

that resonance occurs at radial distance from pulsar B 
of 



1/3 



1/3 



1400 R^Bli^iy, 
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and that the optical depth is given by 



2n ilBTr 



0.2 X{rr)Bli'^,^, 



Here is the radio frequency v 



(8) 

w/(27r) expressed 

in GHz and lobi, — 1.8 x 10^*^5*^12 s^^ is the cyclotron 
frequency at the surface of pulsar B (i.e. for B = B^) 
with 5*12 = 5/(10^^ G). Throughout the paper we set 
the neutron star radius to be 10 km and use i?* as a 
unit of distance. Note that for the A = 1, Tc is approxi- 
mately the ratio of the resonance radius to the radius 
of B's light cylinder, c/^Ib, a result obtained previously 
by Blandford & Scharlemann (1976) and Mikhailovskii 
et al. (1982). 

The determination of the inclination of the binary's 
orbit to the plane of the sky based on the scintillation 
technique (Coles et al. 2004) implies that the radio beam 
of pulsar A passes pulsar B at an impact parameter 
p ^ A X 10^ cm. Since p < r,., at a first glance cy- 
clotron absorption looks like a viable eclipse mechanism, 
although a A ^ 10 — 100 would be required to match the 
eclipse depth. However, a closer look at cyclotron absorp- 
tion, as described below, reveals a problem: A's radio 
radiation heats the particles in B's magnetosphere to rel- 
ativistic energies making synchrotron absorption rather 
than cyclotron absorption the relevant process. 

We denote by F^{ijj,r) the energy fiux per unit fre- 
quency from pulsar A at distance r from pulsar B. A 
nonrelativistic electron or positron absorbs and emits en- 
ergy at rates (Rybicki & Lightman 1979) 

1 
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Here 7 = (1-/3^)^/^ is the Lorentz factor of the electrons 
and positrons which must be close to unity for cyclotron 
absorption and emission to pertain. Balancing by 
E^ yields 

1/2 

(11) 



F^{uB{r),r) 



Lyne et al. (2004) measure a time-averaged flux den- 
sity « 1.6 mJy at h'o = wo/(27r) = 1.4 GHz from pulsar 
A. Taking 600 pc for the distance to J0737-3039 one ob- 
tains Fq = F^{ujq) « 10~® ergs cm~^ at the position 
of B which is se par ated from A by about 9 x 10^° cm. 
Using equation Ullfl. we find that at rr ~ 1250i?*i3*i2 
where 1.4 GHz photons get absorbed the electrons and 
positrons have (3j « 25. This means that the electrons 
and positrons are relativistic and that the cyclotron ap- 
proximation is inapplicable. The timescale for particles 
to become mildly relativistic {E rUeC^) due to cy- 
clotron absorption is 

„2 



heat 



(r) 



nieC 



,2^3/g2 



(12) 



E+{r) 2t:^ FUujB{r),r)' 

which is about 5 s at the position where 1.4 GHz photons 
are resonantly absorbed. 
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A related example of the heating of particles to rel- 
ativistic energies through resonant cyclotron absorption 
of radio waves is given in Lyubarskii & Petrova (1998). 
In their example the pulsar's radio waves heat particles 
streaming along its open field lines. 

3. SYNCHROTRON ABSORPTION. 

The mean cross-section for synchrotron absorption by 
an isotropic distribution of particles with energy jrUeC^ 
is (Rybicki & Lightman 1979) 



34/3r(l/3) S 



5/3 



(13) 



provided wb/i ^ ^ I'^^B- Here T{x) is a complete 
gamma function. 

In the absence of published measurements of pulsar A's 
radio spectrum, we assume that it is a power law with 
most of the energy concentrated at low frequencies (i.e. 
F^^ oc with (5 > 1) as is typical for millisecond pulsars 
(Kuzmin & Losovsky 2001). The position of the low 
frequency cutoff of the spectrum is not very important 
for our problem (although see 

For the cross-section (I13|l and an incident spectrum 
with ^ > 1, particle heating is dominated by the lowest 
frequencies. This also implies that A's radio beam suffers 
a low frequency cutoff which progresses toward higher 
frequency with increasing depth in B's magnetosphere. 
Consequently, the absorbed spectrum of A's radio flux 
at distance r from pulsar B takes the form 



F^(w,r) = Fo j exp[-r^^(w,r)] , 



(14) 



where T^{uj,r) is the frequency-dependent optical depth 
at r. Using equations (j2Jl & H13|) one finds 
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where in arriving at the last expression we have assumed 
that X{r){LUB /l)^^^ is a steeply decreasing function of r 
and 
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dlnr 



(16) 



3 dlnr 

Expression (|15|) is valid provided lu > lob/j] below 
this frequency the relativistic plasma becomes transpar- 
ent. Since extinction at observed frequencies occurs at 
r ~ 10^ cm <C r^^B = c/VLb ~ 10^" cm, a large A(r) is 
required to account for > 1. 

Next we give a simplified evaluation of the energies 
to which magnetospheric particles are heated assuming 
6 = 2 and A independent of r.^ We define a local cut- 
off frequency cji(r) such that (wi (r) , r) = 1. Using 
equation ((T^ and dropping constant coefficients we find 



uji{r) 



7(r) 



A(r 



Hbt 



3/5 



(17) 



Appendix ^contains a more detailed derivation. 



Clearly wi(r) increases as r decreases. Figure illus- 
trates how the low frequency part of Fuj{uj,r) erodes 
with increasing depth in B's magnetosphere.^ The syn- 
chrotron heating rate, E^, is is given by 



Elir) « a{LU,{r))LU,{r)F^{u;,{r),r) 
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which reflects the dominance of radiation with lu ^ 
uJi{r). Balancing the synchrotron heating rate by the 
synchrotron cooling rate H10() , we obtain 



7(r) « 5 X 10^—^ 

-8/5 



A- 



^IbR* 



r 



^0 
52/5 



52/5 



2.5 X 10^ A-«/'^ f'^') B:,\ , 



(19) 



The calculation of 7(r) tacitly assumes that synchrotron 
absorption of A's radio emission is the only source for 
heating particles in B's magnetosphere [cf. Lyutikov & 
Thompson (2005)]. 

To evaluate the size of the eclipsing region re(cj), de- 
fined as a distance from pulsar B at which T^{Lu,re) = 1, 
we combine equations (|17|l and (|19|) to arrive at 



r,{Lu) = 0.54 (^^) 



Fn 



UJQ 



^tJB^N|5/64 

5/64 



^bR* 



meUJB^ XiOB 

380 R^X^'/^^'/'^Blli'' 



11/64 



(20) 



The weak frequency dependence of Tg is an artifact of the 
assumed constancy of A. 

The observed eclipse duration corresponds to ~ 10^ 
cm &l V = 1 GHz. From equation (|20|l we find that 
this implies A « 270. For this value of A, equation H19|l 
yields 7 « 3.2 -B7i2 = 10^ cm so our assumption of 
synchrotron absorption is marginally valid. The optical 
depth cannot be arbitrarily high because radio pho- 
tons of frequency oj are not absorbed within the radius 
at which u ~ ^bH — ^- This "saturation" of t^^ can 
be used to probe the dependence of A on r (see equation 
j32| for a particular example). 

4. ENHANCEMENT OF PARTICLE NUMBER DENSITY IN THE 
MAGNETOSPHERE OF PULSAR B. 

We hypothesize that there is a continuous supply of 
energetic particles within the closed field line region of 
B's magnetosphere from a source located close to the star 
where the magnetic field is strong. Radiation damps the 
gyrational motions of particles born in this region be- 
fore they can slide out along magnetic field lines. Under 
normal circumstances each particle would loop along the 
closed field line on which it was born and strike the star's 
surface in the opposite magnetic hemisphere. However, 

* The low frequency cutoff of A's spectrum is unimportant inside 
the radius at which the cutoff frequency falls below i.i;i(r). 
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illumination by A's radio beam can re-excite the gyra- 
tional motions of particles that stream far away from B 
and then keep them suspended between conjugate mirror 
points (see M 4.2(1 . Thus, both particle heating, making 
synchrotron absorption possible, and particle trapping, 
leading to a higher optical depth eclipse, are mediated 
by A's radio beam. 

4.1. Lifetimes of Energetic Particles. 

After illumination by A's radio beam ceases, parti- 
cles trapped by magnetic mirroring in B's magnetosphere 
cool and are lost. Synchrotron emission reduces their en- 
ergies to TOgC^ on a timescale, tcooi independent of their 
initial 7^1. From equation (|A2|I . we estimate 

W-5|^« 350 s (21) 

As a consequence of relativistic beaming, components 
of momentum parallel and perpendicular to B decay at 
the same rate. In the nonrelativistic regime, velocity 
components perpendicular to B decay exponentially on 
timescale tcooi while those parallel to B remain nearly 
constant. As a result, the trapped particles sediment 
onto the neutron star's surface on timescale tcooi leav- 
ing only the minimum number density of particles, uqj, 
needed to support the corotation electric field. Thus the 
magnetosphere of pulsar B has to refill each time the 
radio beam of pulsar A illuminates it. 

4.2. Trapping of Particles by A 's Beam. 

The particle cooling time, tcooi, given by equation H21|l 
is comparable to the local travel time, ~ r/(/3||c), at the 
distance 

rcooi « l6QR,/3-^^'' Bli'^ (22) 

from pulsar B. Without heating, a particle arriving at 
re from a source located below rcooi would be in its low- 
est gyrational state. To become trapped, it has to gain 
enough perpendicular momentum during its outer mag- 
netosphere passage to mirror above rcooi^ because of adi- 
abatic invariant conservation. 
The minimum gain in required for trapping is 

Api„.„>(/3||7)^(-ec)^(^)', (23) 

where is evaluated at Tg. So the minimum energy, 
^Emim that a particle has to absorb at r ^ rc is 

A P > ^il^ 2 f rcoolY 
2 V '^e / 

«2xlO-X^-eC^^!il(^)". (24) 

Cyclotron absorption of A's unattenuated radio beam 
in a single passage of rc results in an energy increase 

f ■ , s dr 
AE+ ^ / E+{r) — 

^ The expression for the residence time, ~ r/(/3||c), and that for 
rcooi given by equation I22i hold in both relativistic ^ 1) and 
non-relativistic (/Sy <C 1) regimes. 



where the last line is evaluated for S — 2. Use of the 
cyclotron absorption formula ((HJ is appropriate because 
particles streaming past r^ from a source interior to rcooi 
are at most mildly relativistic (see ti44.3() . 

Comparing AEmm with AE^ we find that trapping via 
cyclotron absorption and subsequent mirroring is possi- 
ble for particles with 

Pl'^'l < 0.12 B:T[^y- (26) 

Thus particles reaching r^ ~ 10^ cm with /?|| < 0.4 can 
be trapped. This threshold is not very restrictive (see 
ti44.3|l . so the rate at which A can grow is primarily de- 
termined by the rate at which particles are injected by 
the source near B's surface.^ 

4.3. Source of Particles. 

We can only speculate about possible sources of parti- 
cles in the corotating magnetosphere of pulsar B. The re- 
quirement that B's magnetosphere fill with an appropri- 
ate density of plasma while it is illuminated by A's radio 
beam is not demanding. The total number of particles 
needed to provide the corotating charge density in B's 
magnetosphere is TV^^ {VIbB^,RI/ ec)\TL{rsol R^) ~ 
10"^". Boosting this number by a factor ^ 10^ in 10'^ 
s implies a trapping rate N ~ 10^^ s^^. If every 
trapped particle were born with the energy fimeC^ (sub- 
sequently lost as synchrotron radiation) and only a frac- 
tion /2 < 1 of them were trapped, the source power would 
be P ~ 10^3 /1//2 ergs s-\ For /i = 10 and = 0.1, 
this amounts to ^ 10^^ ergs s~^, much smaller than the 
spin-down luminosity, ~ lO'^" ergs s~^, estimated for pul- 
sar B by Lyne et al. (2004). 

Our favored source is the creation of pairs on closed 
field lines by gamma rays with energies 100 MeV 
emitted by particles accelerated in the outer gap of pul- 
sar B (Cheng, Ho, & Ruderman 1986a,b; Wang et al. 
1998).^ Photons from the outer gap can enter the coro- 
tating magnetosphere and propagate at significant angles 
to magnetic field lines. This facilitates Bj pair creation 
close to the neutron star's surface. 

Pairs born relativistic initially lose energy through syn- 
chrotron radiation on a timescale T'^ - lO'^^ g ^-2 
(see j5]). As a consequence of relativistic beaming, their 
pitch angles remain constant until they become transrel- 
ativistic. Subsequent cooling by gyrosynchrotron radia- 
tion completely damps gyrational motions but preserves 
the component of velocity parallel to the magnetic field. 
Thus particles streaming away from the neutron star 
have /3i| < 1. As we discussed in t|44.2l this allows for 
their efficient trapping when they are illuminated by A's 
radio beam. 

5. STEADY STATE MAGNETOSPHERE. 

Under constant illumination by A's radio beam, B's 
magnetosphere would achieve a steady state. We have 
already described how particle energies would be set by 

^ The two-stream instability might in principle assist in the trap- 
ping of particles, but an unrealistically high number density is re- 
quired for its growth time to be comparable to the particle residence 
time in the magnetosphere. 

^ We are grateful to Jonathan Arons for drawing our attention 
to this possibility. 
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Fig. 1. — Evolution of pulsar A's radio spectrum due to syn- 
chrotron absorption in B's steady state magnetosphere assuming 
X = 2. Solid curves are labeled by distance from B. Dashed curve il- 
lustrate s m aximally absorbed (with optical depth equal to Tmax (t^), 
see eq. 1321 ) spectrum close to B. 



a balance between the absorption and emission of ra- 
diation. Here we concentrate on how the particle en- 
hancement factor A would be fixed in the context of the 
discussion in 21 

The efficiency of particle trapping drops sharply as 
the optical depth at the local cyclotron frequency, 
Tui i^B i'l') , 1^) J increases above unity. The energy that 
fresh particles acquire by cyclotron absorption in pass- 
ing through the magnetosphere is reduced by a factor 
exp[—Tuj{ujB{r), r)] compared to AE+ (see equation 
Thus increasing (w^ (r) , r) lowers the maximum at 
which particle trapping can occur (see H44.2II . The expo- 
nential dependence of trapping on optical depth suggests 
that t^{lob, r) does not greatly exceed unity and also that 
it depends at most logarithmically on r. We set 



(27) 



where x is a parameter we treat as independent of r 
and whose precise value depends upon the rate at which 
trapped particles are lost by unspecified relaxation pro- 
cesses. 

For oj > ujb/j vfG find using equation H15|l that 
Ti^(uj,r) — Ti^{ujB,r){ujB /ujY^^^ . Thus, in the steady state 
magnetosphere photons of frequency uj are absorbed at 
distance re(w) such that ojBire) = lijx~^^^- This gives 



V CO 



1/3 



» 1400 R.x'^^Bliliy-^/'' (28) 

where rr(a;) is defined in equation Q. In steady state 
the eclipse duration would scale as v'^/"^. 

The structure of the steady state magnetosphere for an 
arbitrary power law spectrum of incident radio radiation 



is described by equations IjBlj) and (jB2|) in Appendix IbI 
With typical parameters for J0737-3039 and 6 — 2 we 
obtain 



7(r) w 1.4x' 



-24/55 



r/R. 
103 



36/11 



B 



-12/11 



/ /n \ 49/11 



-20/11 
12 



(29) 
(30) 



Absorption by the extra plasma of A's radio radiation re- 
duces the efficiency of particle heating thus lowering the 
value of 7. The dispersion measure variation during the 
eclipse caused by extra plasma within the steady state 
magnetosphere is at least 2 orders of magnitude below 
the current upper bound of 0.016 pc cm^'^ (Kaspi et al. 
2004). 

The important distinction of the steady state model 
is that it predicts A. This determines the lower and up- 
per frequencies between which synchrotron absorption is 
effective at a given r: 

/ /D \ -69/11 

-1 « 2 GHz v^^/^M ^ ) b'^' 



2tt 



103 ; "-^12 . 
/ X 39/11 

« 5.6 GHz X--/- (^) B:,T' (31) 

It follows that the maximum optical depth at frequency 
ly, an observable quantity, is given by 



,(z.)«3.2x^^^^^i^, 



15/23,-20/23 



(32) 



Figure depicts the evolution of the spectrum of A's 
radio beam with depth in B's magnetosphere. Power 
at low frequencies is gradually eaten out by synchrotron 
absorption as the beam propagates deeper into B's mag- 
netosphere. 

Illumination of B's magnetosphere by A's radio beam 
probably starts only short time prior to eclipse. So 
we may be witnessing radio beam attenuation by dy- 
namically evolving plasma in B's magnetosphere. The 
timescale for cold particles to become trans-relativistic 
via cyclotron absorption evaluated from equation H12|l is 
rather short, typically < 10 s. Energies of relativistic 
particles rise exponentially via synchrotron absorption 
on timescale 



t 



heat 



0.78 



:F^{ujb) 



320 s 



r/R, 

103 



(33) 



so long as the radio spectrum of pulsar A remains unat- 
tenuated. Heating has an exponential character for 5 = 2 
because as 7 grows particles can absorb the incoming 
radio photons at lower frequencies because ujbH de- 
creases making more energy from the unabsorbed spec- 
trum available to heat them. This partly compensates 
for the decrease of the absorption cross section with in- 
crease of 7 (see eq. QSl)- Once 7 grows so large that 
TuiiuJB/j,r) becomes comparable to unity, heating is less 
efficient and theat oc 7^/3 estimate that particles 

reach 7 10^ at r « 10^ cm on a timescale of ~ 5 min. 
This is comparable to estimates of several tens of min- 
utes for the duration of the illumination period judged 
from the orientation of pulsar A's spin and dipole axes 
suggested by Demorest et al. (2004). Not much can be 
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said about the timescale needed for the particle density 
to reach its steady state value. Presumably this is largely 
controlled by the rate at which the source near pulsar B 
is able to supply fresh trans-relativistic particles since 
these are readily trapped. 

6. REPROCESSED RADIATION. 

Energy absorbed by particles in B's magnetosphere 
from A's radio beam is reemitted as synchrotron radi- 
ation, albeit with some time delay. Absorption takes 
place at w ~ wi, see equation H17(l . close to the mini- 
mum frequency possible, lobIIi and is reemitted at the 
considerably higher frequency u ~ 7^073. This reemitted 
radiation, although unimportant locally, may dominate 
the primary radiation from pulsar A deeper in pulsar B's 
magnetosphere. 

Here we attempt to calculate the properties of the re- 
processed radiation for the steady state magnetosphere 
by applying results from SjSl This requires choosing a 
low frequency cutoff, Vmim for the assumed power law 
spectrum of A's radio emission. Because there is little 
evidence for low-frequency cutoffs in the spectra of mil- 
lisecond pulsars above 100 MHz (Kuzmin & Losovsky 
2001), wc normalize Vmin to this frequency. From equa- 
tion (|28(l we deduce that the energy carried by pho- 
tons with V ~ v„iin is absorbed at a distance r^ax ~ 
— 1/3 

3500i?^i^^j„ g from pulsar B and reemitted at frequency 

l^max ~ l'^{rmax)^B{rmax)l'i'^ ~ GIHz. The 

latter represents the upper cutoff frequency of the repro- 
cessed radiation because ^'^ojb increases with increasing 
r (see UHl)- A detailed calculation shows that radiation 
reemitted by all the relativistic particles within rmax has 
a flat power law spectrum F„ with index close to zero. 
The local flux of reprocessed photons at Umax can be 
estimated from F^[vmax)vmax ~ F^{.Vrmn)Vmin yielding 
Fv{vmax) ^ 2 X 10~^ ergs cm~^ s~^ Hz"-'^ for a primary 
spectrum with 6 — 2 and time-averaged intensity (at 
B's location) F^{IA GHz) = 6 x 10"^ ergs cm'^ s'^ 
Hz~^. Reprocessed spectrum intersects the primary ra- 
dio spectrum of pulsar A at ~ 50 GHz and at frequencies 
higher than that one has to solve the radiation transfer 
problem to determine particle heating. Luckily, this is 
far enough from 1 GHz region of spectrum where eclipse 
observations are usually taken so that we do not need 
to worry about such complications. Unfortunately, the 
reprocessed radiation can hardly be detected on Earth 
because of very small covering fraction of B's magneto- 
sphere as seen from A. 

7. DISCUSSION. 

A simple dynamical picture of pulsar A's eclipse by 
pulsar B's magnetosphere emerges from our considera- 
tions. Prior to the onset of illumination by A's radio 
beam, the number density in the corotating part of B's 
magnetosphere is equal to ncj- When A's radio beam 
strikes the magnetosphere of pulsar B, perhaps ^ 10 
min prior to B's inferior conjunction, particles initially 
present in the magnetosphere rapidly (in < 10 s) become 
relativistic. Their energies continue to rise until either 
synchrotron absorption is balanced by synchrotron emis- 
sion or illumination by A's radio beam ceases. At the 
same time neutral plasma accumulates in the magne- 
tosphere as the result of the trapping of particles that 



stream out in cold beams from a source near B. It is 
unclear whether the magnetosphere reaches steady state 
prior to the eclipse or whether is still evolving. In ei- 
ther case a density enhancement A ^ 200 — 300 would 
be required to match the observed depth and duration of 
the eclipse. Some time, perhaps several tens of minutes, 
after the eclipse illumination of B's magnetosphere by 
A's radio beam ceases and the particles cool on a typical 
timescale tcooi ^ 10^ — 10^ s given by equation H21|l . As 
a result particles no longer mirror and all plasma beyond 
that needed to maintain the corotation electric field is 
lost. This entire cycle repeats each orbital period. 

We speculate that part of the asymmetry between 
eclipse ingress and egress reflects a rise in plasma density 
in B's magnetosphere during the 30 s the eclipse lasts. 
The optical depth, is very sensitive to increasing A. 
Not only is directly proportional to the total column 
density of absorbing particles, it is also proportional to 
the absorption cross section per particle which varies as 
^-5/3 j^j^^ ^ decreases with increasing A. Thus if A were 
growing on a timescale of minutes during eclipse, egress 
would be deeper than ingress and the eclipse centroid 
would occur slightly after B reached inferior conjunction 
since Te grows as A increases (see eq. (201). In this picture 
the smooth ingress may partly be caused by the time- 
variable optical depth at a fixed location rather than the 
variation as the radio beam samples smaller values of 
r. On the contrary, during egress the optical depth is 
higher and the abrupt termination of eclipse may signal 
the emergence of the radio beam from behind an almost 
opaque screen. This explanation requires proper timing. 

Our model predicts a variation of the polarization of 
A's radio emission during the course of the eclipse. The 
polarization signal should be strongly correlated with 
B's rotational phase in a manner similar to the eclipse 
lightcurve variations found by McLaughlin et al. (2004b). 
This prediction stems from the fact that synchrotron ab- 
sorption of photons in different polarization states is sen- 
sitive to the angle between the photon k vector and the 
direction of the magnetic field. Modeling the polariza- 
tion signature would be facilitated because B's magnetic 
field should be nearly dipolar at since i?* ^ <C r^o- 

Effects of general relativity are very important in 
J0737-3039. Lai & Rafikov (2004) have demonstrated 
that gravitational light bending can significantly (by 
~ 30%) change the minimum impact parameter at which 
A's radio beam passes B. Thus gravitational lensing 
must play a significant role in shaping the eclipse pro- 
file. Moreover, because of the binary's finite orbital ec- 
centricity, e ~ 0.088, periastron precession driven mainly 
by the effects of general relativity forces a 21 yr periodic 
variation, from a \ cos i | ( 1 — e) to a | cos i | ( 1 — e) neglecting 
lensing effects, of the pulsars' minimum projected sepa- 
ration on the plane of the sky (Burgay et al. 2003). By 
changing the minimum impact parameter of A's radio 
beam with respect to B, this should produce observable 
eclipse profile variations. 

Lyutikov (2004) and Arons et al. (2004) have sug- 
gested that the eclipse of A's radio beam is caused by 
synchrotron absorption in the magnetosheath that forms 
when A's relativistic wind impacts B's magnetosphere. 
Detailed polarization observations offer a means to dis- 
tinguish this viable alternative from our model. 
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The model we have presented is necessarily rather sim- 
plistic — it is one-dimensional, it assumes that particle 
distribution functions are isotropic, and so on. Where 
possible, future work should relax these constraints. A 
quantitative description of the particle source giving rise 
to the plasma density enhancement might also be pur- 
sued. Further observations have the opportunity to re- 
veal additional clues to properties in the eclipse region. 
Knowledge of A's radio spectrum is crucial for calculat- 
ing particle energies, trapping efficiencies, and the eclipse 
duration. Time-resolved polarization of A's radio emis- 
sion during eclipse would constrain the magnetic field 
geometry and particle distribution anisotropy. 

Clarifying details of this remarkable example of nonlin- 
ear coupling of relativistic plasma to external radiation 



may provide clues for understanding other puzzling phe- 
nomena such as the modulation of B's radio emission by 
radiation from A (McLaughlin et al. 2004a) and the peri- 
odic variations of the pulsar B's brightness that correlate 
with its orbital phase (Ransom et al. 2004). 
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APPENDIX 

PARTICLE ENERGY. 

To properly compute the heating of particles by synchrotron absorption we evaluate 



1 



Kir)-!, 



F^{uj,r)(Ts {ijj)duj 



4^^ 



j/7 



34/3r(l/3) B 



(Al) 



where the factor 1/2 roughly accounts for shadowing as B rotates. The lower limit of integration is extended to 
zero (instead of ljbIi) because we are assuming that T^^itos /l ^ 1. The largest contribution to the local heating 
comes from lo such that T^(w,r) ^ 1 which, coupled to the condition that lob/i < ^ and equation (|15|) . assures that 
Tu{i^B/l,r) > 1. 

Substituting from equation H15() into equation IjAip and calculating the integral over duo we find 

n -(2+35)/5 



An 



34/3Cir(l/3) 



^^'^ ^^^^ CUJbFq f UJq'^ ^ ^ 



5 V 5 

Balancing this heating rate by the cooling rate (|10|l yields 
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Substituting equation ljA3|l into equation H15|) . assimring t^{u!) — 1, and solving for r^, we obtain the size of eclipsing 
region at frequency lu: 



re(w) = Ri, 



47r 



34/3r(l/3)Ci 



3/32 



20 



277rCir(2i^) 



5/64 ,(3_„ 



^^0 



^0 
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-5/64 



11/64 



(A4) 



The constant coefficients in these expressions depend upon the values of the parameters S and Ci- For the radio 
spec trum of pulsar A with 5 = 2 and A independent of r, it follows that Ci = 64/3. Then the coefficient in equation 
(|X3)l is « 481. 

STEADILY ILLUMINATED MAGNETOSPHERE. 

Combining equations (|15|l . (|27|l . and (jA3p . 7 is self-consistently determined to be 



7(r) 



97r2 r 



'2+35^ 
\ 5 J 



2 + 3S Ft 



31/35 r(i/3) 

From equations l(T^ . and IjBip we find that 



X 







3/11 
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34/35r(l/3) 



47r 
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3{3-S) 



Fo 



Wo 
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5/11 



(Bl) 



(B2) 



Equations HB1|I and (|B2(I imply that (i = 5, see equation (|16|) . These formulae determine the structure of the steady 
state magnetosphere for arbitrary S > 1. 
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